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We present a new asymptotically-flat time-machine model made solely of vacuum and dust. The 
spacetime evolves from a regular spacelike initial hypersurface S and subsequently develops closed 
timelike curves. The initial hypersurface S is asymptotically flat and topologically trivial. The 
chronology violation occurs in a compact manner; namely the first closed causal curves form at the 
boundary of the future domain of dependence of a compact region in S (the core). This central core 
is empty, and so is the external asymptotically flat region. The intermediate region surrounding the 
core (the envelope) is made of dust with positive energy density. This model trivially satisfies the 
weak, dominant, and strong energy conditions. Furthermore it is governed by a well-defined system 
of field equations which possesses a well-posed initial-value problem. 



I. INTRODUCTION 

This paper deals with the possibility of formation of 
closed causal curves (CCCs) in spacetime, within the 
framework of General Relativity. By causal curves we 
mean either timelike or null curves. The main ques- 
tion at the background is the following: Is it possi- 
ble that CCCs will spontaneously evolve from rather 
"normal", non-pathological, initial conditions? In most 
of the classic solutions of the Einstein equations (e.g. 
Minkowski, Schwarzschild, Robertson- Walker) no such 
CCCs occur, although several examples of spacetimes 
which do admit closed timelike curves (CTCs) are known 
00 II II HUB Hi- We dub such spacetimes, which ad- 
mit CTCs, as "time-machine (TM) spacetimes". How- 
ever, most of these examples suffer from pathologies or 
problematic ingredients which question their relevance to 
physical reality. Following is a list of three basic require- 
ments that eliminate most of the TM models proposed 
so far: (i) The spacetime should admit a regular space- 
like initial hypersurface (a partial Cauchy surface) S; (ii) 
asymptotic flatness; (hi) the weak energy condition Q. 
Thus, Godel's rotating-dust cosmological model vi- 
olates conditions (i,ii), as does Tipler's rotating-string 
solution 01; the wormhole model by Moris, Thorne and 
Yurtsever violates condition (iii). Gott's solution 
of two infinitely long cosmic strings violates condition (ii) 
1C|, and Mallett's solution violates condition (i) (see 

A few of the previous models 000 do satisfy the 
above requirements (i— iii) . These are all asymptotically 
flat, topologically trivial models, in which the CCCs are 
born inside a certain torus. These models, however, fail 
to satisfy some other basic requirements discussed below. 
In particular, the models 00 violate the strong energy 
condition [9( ■ Here we present a new model which better 
satisfies all these requirements.. 

In considering the physical relevance of a TM model 
(like any other kind of spacetime model), one of the most 
important aspects is the physical suitability of the space- 
time's matter content. One certainly would like to im- 



pose the weak energy condition, and preferably also the 
strong and dominant energy conditions. But com- 
pliance with the energy conditions is not enough. One 
would also like the spacetime's energy-momentum tensor 
to coincide with some known matter field, governed by 
a well-known, well-posed field equation, for at least two 
obvious reasons. First, if we cannot associate the energy- 
momentum tensor with some known matter field, then we 
cannot tell whether this energy-momentum distribution 
can be obtained in reality. Second, the following question 
concerning spacetime dynamics is at issue: Can one de- 
sign a "normal" initial configuration such that the laws 
of evolution will subsequently force the spacetime to de- 
velop CCCs and to violate chronology? A system (space- 
time -(-matter) which does not admit a well-defined set 
of evolution equations will be inadequate for addressing 
such a question. 

Motivated primarily by the last argument, in Ref. 
we presented a TM model in which the initial hypersur- 
face is composed of three parts: An external asymptoti- 
cally flat region, the "envelope" (an intermediate region), 
and a compact toroidal region at the center, to which we 
refer as the "compact core". In this model CCCs evolve 
inside the compact core, in a manner which is causally in- 
dependent of the surrounding regions. The central com- 
pact core, and also the external asymptotically-flat re- 
gion, are vacuum, but the envelope is made of matter. 
This matter satisfies all three energy conditions men- 
tioned above, yet the envelope's energy-momentum has 
not been recognized as any known form of matter field. 
Nevertheless, since the internal core is made of pure vac- 
uum, and the formations of CCCs inside the core is guar- 
anteed independently of the envelope's evolution, this 
model fulfills its main goal at least to some extent — 
it successfully demonstrates how the laws of spacetime 
dynamics inevitably lead, in a certain situation, to the 
violation of chronology — provided that the initial config- 
uration of energy-momentum at the envelope could be 
realized by some real matter field. But still the question 
remains whether such a matter field exists or not. One of 
the main goals of this paper is to address this difficulty. 
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The spacetime model constructed here is similarly 
composed of three parts: An internal vacuum core, 
an external asymptotically flat vacuum region (the 
Schwarzschild geometry), and a non-empty intermedi- 
ate region (the "envelope"). Here, however, the enve- 
lope's matter will be simply dust (namely, a perfect fluid 
with zero pressure), with non-negative energy density. 
This kind of matter trivially satisfies the weak, dominant, 
and strong energy conditions, ft also yields a well-posed 
initial- value problem |9J. Dust is probably not the most 
realistic or fundamental description of matter, yet it has 
been proven useful in addressing various issues of princi- 
ple in General Relativity, e.g. the dust Robertson- Walker 
cosmology, the Oppenheimer-Snyder model [13] of homo- 
geneous dust collapse (the first model to demonstrate the 
formation of a black hole in gravitational collapse), and 
the formation of naked singularities in spherical dust col- 
lapse [131 ]. In the last two problems, significant progress 
was first made by exploring a dust model, but qualita- 
tively similar results were subsequently observed in mod- 
els with nonvanishing pressure (see e.g. [H[). In fact, 
it appears that in our problem as well it will not be dif- 
ficult to generalize the present dust model to a perfect 
fluid with pressure, but this is beyond the scope of this 
paper. 

The present model also differs from that of Ref. [7[ 
in the type of vacuum metric employed for the com- 
pact core. In Ref. Q we used a vacuum solution lo- 
cally isometric to a pp-wave [HI ] spacetime. Here we 
use a vacuum solution locally isometric to a "pseudo- 
Schwarzschild" geometry (namely, one obtained from the 
Schwarzschild geometry by a Wick rotation), which we 
describe in Sec. IIII1 One of the differences is that the 
pseudo-Schwarzschild core metric can be easily repre- 
sented in a diagonal form, which globally covers the core 
metric from the initial hypersurface and up to the Cauchy 
Horizon (CH). We are not aware of such a global diagonal 
representation of the core metric of Ref. [3] • Another dif- 
ference is that the present core metric admits an initial 
hypersurface with an especially simple extrinsic curva- 
ture, as described in section Hill — again, we are not aware 
of such a possible choice of initial hypersurface in the core 
metric of Ref. 0]. Both factors, the diagonal core metric 
and the simple form of the extrinsic curvature, greatly 
simplify the construction of the initial data for the TM 
model with dust envelope. 

There is another, more significant, difference between 
the two core metrics: The locally pseudo-Schwarzschild 
metric is similar to the (four-dimensional version of the) 
Misner space [16], as its CH is entirely generated by 
closed null geodesies (CNGs). In the locally-pp metric, 
on the other hand, the CNGs are generically isolated. 
The pseudo-Schwarzschild core metric is also similar to 
the Misner space in that two non-equivalent analytic ex- 
tensions beyond the CH exist [TtJ ■ On the contrary, only 
one possible extension beyond the CH appears to exist 
in the locally-pp core metric of Ref. Q • 

The Misner-like form of our present core metric has 



both advantages and disadvantages. A CH generated by 
CNGs is claimed [HI to be unstable against "fragmen- 
tation" into a set of isolated null geodesies. This is a 
disadvantage which might motivate us to try construct a 
dust envelope for the locally-pp core of Ref. Q as well, 
but this is beyond the scope of the present paper. At any 
rate we do not attempt to address issues of stability in 
this paper. 

But the Misner-like core used here also has advantages. 
In our previous model [7[, the question arises whether 
the closed null geodesic N at the CH is adequately pro- 
tected against a singularity which might form at the fu- 
ture boundary of D + (S) and approach arbitrarily close 
to N [19]. Originally we thought that the CNG N is 
causally protected against this scenario, for the following 
reason: The initial hypersurface S has a compact core 
So, such that the CNG N is located at the boundary of 
D_|_(So). The vacuum solution in the entire region D + (So) 
is known explicitly, and is regular throughout. This en- 
sures that no singularity can form at (the boundary of) 
D + (So) and endanger the regularity of the neighborhood 
of N. However, the domain D + (S) is larger than D+(So)- 
Although the boundaries of these two domains coincide 
at N, it may well be the case that there is a separation be- 
tween these two boundaries, and the set D + (S)— D + (So) 
gets arbitrarily close to N. This structure may allow the 
possibility of a singularity which evolves at the future 
boundary of D + (S) (but outside D + (So)) and extends 
arbitrarily close to N. In that case, N would be a regu- 
lar CNG (with a regular neighborhood) as viewed from 
D + (So), but would still lack a regular neighborhood in 
D+(S). This might be harmful for an extended (test) ob- 
server attempting to cross the CH (of S) through one of 
the points on N, in order to penetrate into the region of 
CTCs. 

It is not clear if this potential problem is realized in the 
model Q. It is hard to say, because the exact solution 
for the time-evolving metric is only known in the inter- 
nal vacuum core and in the external region, not in the 
envelope. Therefore, we do not know the full structure of 
D+(S), and, in particular, what kinds of singularities it 
develops, if any. Fortunately this potential problem does 
not apply to the present model. As will be shown in Sec. 
IIIIl due to the Misner-like form of the CH, the bound- 
aries of D + (S) and D + (So) do overlap in a set denoted Hi 
below. This set, which is a portion of the CH, includes a 
continuum of CNGs. A sufficiently small neighborhood 
of any such CNG, restricted to D + (S), is entirely con- 
tained in D + (So). Since the metric throughout D + (So) 
is known explicitly and is perfectly smooth, no singular- 
ity which might evolve at (the boundary of) D + (S) can 
get close to any of these CNGs. 

As previously stated, the main underlying question is 
the possibility of triggering the onset of chronology vi- 
olation. In other words, is it possible to design initial 
conditions for which the laws of dynamics will inevitably 
lead to violation of chronology? But there is a build-in 
logical difficulty in the formulation of this question: If 
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indeed CTCs form, then the portion of spacetime con- 
taining the CTCs is by definition outside the future do- 
main of dependence of any initial hypersurface S. In what 
sense can one then state that the chronology violation 
has "emerged from the initial conditions on S"? This is 
indeed a difficulty, but nevertheless we propose a set of 
conditions which, when satisfied, provide meaning to the 
statement that the violation of chronology was triggered 
by the initial conditions on S. These conditions are: (i) 
H + (S) contains CNGs (therefore the Cauchy evolution of 
the initial data on S unambiguously leads to some sort of 
chronology violation); (ii) The analytic extension of the 
metric beyond H + (S) [l7| (or some portion of H+(S)) in- 
cludes CTCs in the immediate neighborhood of H+(S); 
(iii) Any smooth extension of the metric beyond H+(S) 
(or some portion of H+(S)) will include CTCs in the im- 
mediate neighborhood of H + (S). 

As a simple application of these criteria, consider the 
analytically-extended geometry of a Kerr black hole. 
This spacetime is known to admit CTCs deep inside 
the black hole. The CTCs are located beyond the inner 
horizon — a null hypersurface which serves as the CH for 
any initial hypersurface in the external universe. We shall 
not regard this spacetime as a " time-machine model" as 
it fails to satisfy any of the above criteria (i— iii) . In par- 
ticular, the inner horizon does not contain CNGs. On 
the other hand, the model presented in this paper does 
satisfy all three conditions. 

When addressing the possibility of constructing a time 
machine, one would primarily be interested in the situa- 
tion where the construction process takes place in a finite 
region of space. A simple criterion which captures this 
idea is the following: We shall say that the time machine 
is compactly constructed if the initial hypersurface S in- 
cludes a compact set So such that the Cauchy evolution 
of the initial data on So leads to chronology violations; 
That is, the closure of D + (So) includes CCCs (specifically 
this means that H + (So) includes CNGs). 

Hawking [l8| earlier introduced a different notion of 
compactness called compact generation. A CH is said to 
be compactly generated if all its null generators, when 
past-propagated, enter a compact region of spacetime 
and never get out of it. This criterion differs from 
the notion of compact construction formulated above. 
The time-machine model presented in this paper, as well 
as our previous models QHuIj are all compactly con- 
structed but might not be compactly generated. 

The above discussion led to several criteria which one 
might apply to any candidate model attempting to de- 
scribe the process of "constructing a time-machine" in 
our physical spacetime. In the next subsection we collect 
these criteria and list them in a more systematic man- 
ner. Then section HT1 outlines the structure of our model 
spacetime and its various parts (central core, envelope, 
and external region). In Sec. IIIII we present our core 
metric and discuss its main properties. Section IIVI out- 
lines the initial- value setup, and the constraint equations 
which must be satisfied by the initial data. Then Sec- 



tions [V] and I VII describe the construction of the desired 
initial data (3-metric and extrinsic curvature) on the en- 
velope and external parts of the initial hypersurface S, 
respectively. In Sec. IVIII we summarize and discuss some 
of the problems and open questions remaining for future 
research. 



A. Criteria for a physical time-machine model 

Here we collect the various criteria which emerged in 
the discussion above (plus one more criterion related to 
the space topology). It should be emphasized that we 
do not attempt here to postulate a strict, formal defini- 
tion of a "TM model". Rather, our goal here is to list 
the various criteria which we find relevant. This list may 
serve as a useful basis for discussing the physical rele- 
vance of various models which attempt to describe "TM 
construction" : 

1. The spacetime should admit a spacelike initial hy- 
persurface (a partial Cauchy surface) S; 

2. The initial data on S should be sufficiently regular; 
Namely, both the spatial 3-metric and the extrinsic 
curvature should be for a sufficiently large k. 
(k should be, say, 4 or larger in order to guaran- 
tee a well-defined time evolution. The construction 
below yields C'- 00 - 1 initial data.) 

3. Asymptotic flatness, 

4. The spacetime's matter content should satisfy the 
Energy conditions. This may be divided into two 
categories: (4a) The weak energy condition, and 
(4b) The dominant and strong energy conditions. 

5. The causal connection between S and the chronol- 
ogy violation: (5a) H + (S) should contain CNGs; 
(5b) The analytic extension of the metric beyond 
(some portion of) H+(S) should include CTCs in 
the immediate neighborhood of H + (S); (5c) Any 
smooth ("hole-free" [H|@|) extension of the met- 
ric beyond (some portion of) H+(S) should include 
CTCs in the immediate neighborhood of H + (S). 

6. "Causal protection" of the CNG: H+(S) includes a 
CNG N admitting a neighborhood in D+(S) which 
is perfectly regular. 

7. Compact construction: S should include a compact 
set So such that H + (So) contains CNGs. [Further- 
more, the criteria 5,6 above should apply to a por- 
tion of H+(S) which is also contained in H+(So).] 

8. The topology of S should be trivial (R 3 ). 

9. The energy-momentum tensor will correspond to 
a known matter field, which yields a well-posed 
initial-value problem. This is especially crucial for 
the core metric inside D + (Sq) (which itself develops 
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chronology violation), but is also desired (though 
perhaps to a lesser extent) for the outer parts of 
the time- machine model. 

We may also add the following, wider (and loosely for- 
mulated) requirement concerning the spacetime matter 
content: 

10. The matter field will be as elegant and/or realistic 
as possible. 

From the point of view of classical General Relativity, the 
most elementary and elegant type of energy-momentum 
tensor is obviously the vacuum, T a p = 0. 

The present model satisfies all these requirements. Our 
previous model [3] does not fully satisfy criterion 9, be- 
cause the envelope is made of an unrecognized matter 
field (though the core is vacuum). Also it is not clear if 
criterion 6 is satisfied by it. Our earlier models also 
fail to satisfy criterion 4b. In fact, none of the previous 
models satisfy the combination of criteria 1,2,3,9. 

II. AN OVERVIEW OF THE SPACETIME'S 
STRUCTURE 

Our model spacetime is composed of three parts: the 
central active vacuum core, the dust envelope, and the 
external asymptotically-flat vacuum region. Correspond- 
ingly, the initial hypersurfacc S (a partial Cauchy surface) 
will be composed of three parts: 

1. The internal vacuum core, located inside a certain 
torus To. 

2. The envelope, an intermediate region located be- 
tween the torus T and a two-sphere R s surround- 
ing it. (R s is the two-sphere which is later denoted 
R=R 2 .) 

3. The external vacuum region, located outside the 
two-sphere R s . (This external region actually cor- 
responds to a certain spacelike hypersurface in the 
Schwarzschild geometry.) 

These three parts of S will be denoted So, Si, and S2, 
respectively. 

In a similar manner we divide the evolving 4- 
dimensional spacetime into three regions. To be more 
precise, it is the predictable portion of spacetime, namely 
D + (S) and its closure, which we divide and associate with 
the various parts of S. The internal region Mo is D + (So); 
the external region M2 is D + (S2), and the intermediate 
region Mi is the intersection of D + (S) with J + (Si). Both 
Mo and M2 are pure vacuum regions: Mo is a compact re- 
gion which includes CCCs and hence constitutes the TM 
core, and M2 is a portion of the Schwarzschild geometry 
(the external part, which extends to spacelike infinity). 
The intermediate region Mi is made of dust (though it 
also includes vacuum parts). 



In both vacuum regions Mo and M2 the 4-geometry is 
known to us explicitly: It is the metric (|1I2[) in region 
Mo, and the Schwarzschild geometry in region M2. In 
the dust region Mi, on the other hand, the evolving 4- 
geometry is not known explicitly. Instead, it is described 
in terms of the corresponding initial data on S (section 
E] below). 

III. INTERNAL CORE METRIC 

We start our construction from the vacuum metric 

ds 2 = -(l-2u/r)- 1 dr 2 + (l-2n/r)dt 2 +r 2 (d9 2 +smh 2 9dip 2 ) , 

(1) ' 

where \i is an arbitrary positive constant. Here 9 takes all 
positive values, whereas (p admits the usual periodicity 
< tp < 2n. This metric is obtained from the standard 
Schwarzschild metric by a Wick rotation 9 — > i9, and we 
shall refer to it as the pseudo-Schwarzschild metric. The 
coordinate t is assumed here to be periodic, 

< t < I, 

where I is a free parameter which we take to be greater 
than some minimal value Z m i n [specified in Eq. (|23|) be- 
low]. 

The above metric has a coordinate singularity at 
r = 2/i, analogous to the Schwarzschild's horizon (later 
we shall remove this singularity by transforming to 
Eddington-like coordinates). At this stage, and through- 
out most of this paper, we shall primarily be interested 
in the range r > 2/i. Note that in this range t is a spatial 
coordinate (i.e. gu > 0) and r is a time coordinate, which 
we take to be past-directed (namely, r decreases on mov- 
ing from the past to the future). Thus, all hypersurfaces 
r = const > 2/i are spacelike. 

To overcome the coordinate singularity at r = 2/i we 
now transform to "Eddington-like" coordinates in the 
usual manner: We define 

v = -(t + r*), 

where 

r* = r + 2/iln(r/2/i - 1) . 

In the (r, v) coordinates the metric becomes 

ds 2 = (l-2u/r)dv 2 + 2dvdr + r 2 (d9 2 +sinh 2 9d<p 2 ) . (2) 

The coordinate v has the same periodicity as that of t, 
namely, for given (r, 9, (p) the point v = I is identified 
with v = 0. Q 

A. Formation of CTCs 

The closed orbits of constant r, 9, cp admit the one- 
dimensional line element 

ds 2 = (1 - 2n/r)dv 2 . 
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These orbits are spacelike throughout r > 2/i, but they 
become timelike — namely CTCs — at r < 2/i. Note that 
the metric ((2]) passes r = 2/i in a perfectly regular man- 
ner, as all components are regular and 

det(g) — —r 4 sinh 2 8 

is nonvanishing there. 

The " hypersurface" r = is a true, timelike, curva- 
ture singularity. Our analysis throughout this paper is 
restricted to the range r > 0. 

The hypersurface r = 2/i is null, and its generators are 
the curves of constant 8, tp, which are all CNGs. This 
hypersurface is in fact the CH for any partial Cauchy 
surface r = const > 2/i. It also serves as the chronology 
horizon for the metric Namely, all points at < r < 
2/i sit on CTCs (e.g. the curves of constant r,9,tp), but 
none of the points at r > 2/i do (because the region r > 
2/i is foliated by the spacelike hypersurfaces r — const). 



namely tp = (2ir/l)t. Then the inner-core 3-metric be- 
comes 

ds 2 = rl(d8 2 + sinh 2 9dtp 2 ) + L 2 d0 2 , (6) 

where 

L = (Z/27r)Vl-2 M /ro = ^-L= . 

Note that the 3-metric (|BJ| is cylindrically-symmetric, 
with 8 serving as the "radial" coordinate and tp as the 
azimuthal coordinate. [Later we embed this 3-metric as 
the core of a global asymptotically-flat hypersurface (sec- 
tions IVIVII below) . This global hypersurface is axially- 
symmetric. It should be clarified that it is <p, not tp, 
which becomes the global azimuthal coordinate.] 

C. Truncating the internal core metric 



B. Initial hypersurface for the internal core 

When discussing the initial-value problem for the 
above spacetime, we shall consider an initial hypersur- 
face located at r > 2/i. It will be convenient to express 
the metric in the diagonal form {1} (the coordinate sin- 
gularity at r — 2/i will not pose any difficulty, as it takes 
place away from the initial hypersurface). 

For any spacelike hypersurface r = const > 2\i, the 
spatial 3-metric is 

ds 2 = (1 - 2u/r)dt 2 + r 2 (d8 2 + sinh 2 ddp 2 ) = h ab dx a dx b . 

(3) 

Hereafter the indices a, b run over the three spatial coor- 
dinates. The extrinsic curvature K\ of such a hypersur- 
face is 



Jf| = 4(l-2 M /r)-V». 



tf£ = r- V(l - > 



and all other components vanish. The two distinct eigen- 
values coincide at r = 3/x in which K h a = kgS^, or 



Kab = k h ab 



where 



27fi 



(4) 



(5) 



The form ^ greatly simplifies the constraint equations, 
therefore we shall take our initial hypersurface (for the 
internal core) to be at 

r = 3(i = r . 

We denote this hypersurface by S. 

For later convenience we transform the periodic coordi- 
nate t of the internal vacuum core into a new coordinate 
(j) with a standard periodicity 

< cj> < 2tt , 



The hypersurface r = 3/i with the spatial 3-metric © 
is not asymptotically-flat. In order to match it to an 
asymptotically-flat exterior, we have to truncate the in- 
ternal 3-mctric at a certain two-surface. Any two-surface 
8 = const > on the three-surface r = is a torus 
(parametrized by the two periodic coordinates tp, tp). We 
shall truncate the initial 3-metric (JBJ on such a 2-surface 
8 = const = do. We denote the portion 8 < 8a of the 
three-surface r = r by So 

Consider now the set D + (So), namely the future do- 
main of dependence of So- This set has a limited extent 
in the time r (because none of the points at r < 2/i belong 
to this set). At 2/i < r < ro, D + (So) will be bounded 
by the null geodesies of constant tp, tp which emanate at 
r = ro from 9 = 9q and propagate towards smaller 9 val- 
ues. [Note that geodesies with different tp, <p will have the 
same orbit 9(r), due to the cylindrical symmetry..] For 
our time-machine construction it is crucial that D + (So) 
will include a portion of the chronology horizon at r = 2/i. 
This demand will impose a minimal value for 8q, as we 
now discuss. 

The above null orbits of constant tp, <p which bound 
D + (Sq) satisfy the differential equation 



d9 
dr 



V~9rr 



1 



\/9ee Ty/1 — 2/i/r 
The general solution of this equation is 



0(r)=21n 07 (2/i) + - 1 



(7) 



where Cg is an integration constant. This constant is de- 
termined from the initial value 8 — 8q at r = 3/x, namely 



where 



C e = 8 



= In 2 + V3 = 1.317. 
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At r — 2/i the term in squared brackets in Eq. ([7]) 
vanishes. We therefore demand Cg > 0. Thus we shall 
take the cutoff value 8q to be > 8 C . This ensures that the 
portion 

< < 9 - 6 C 

of the chronology horizon at r = 2/j will be included 
in the boundary of D + (So). In particular this portion 
includes an open set of CNGs. 



D. The relation between H+(E), H+(So), and 
H+(S) 

In the non-truncated core metric © the CH associated 
with the (complete) initial hypersurface E, H + (E), is the 
hypersurface r = 2fi. Once the core metric is truncated 
(as described above), the structure of the CH changes. 
The truncated part So has its own CH, denoted H+(So). 
In addition the global asymptotically-flat initial hyper- 
surface S (which contains So as its core) has its own CH, 
denoted H+(S). Here we shall briefly discuss the relation 
between these various CHs. 

Consider first the structure of H + (So). From the dis- 
cussion in the previous subsection it follows that this 
hypersurface is composed of two parts: (i) the por- 
tion < 8 < 9q — 8 C of the hypersurface r = 2/.1, 
which we denote Hi , and (ii) a null hypersurface denoted 
H2, associated with the truncation of the core metric, 
whose generators emanate from the truncation 2-surfacc 
r = 3/i, 8 = 8 . These generators follow the orbit (JT]) (for 
each ip and </>). The part Hi is a portion of H+(E). This 
part (unlike H2) is entirely generated by CNGs. 

The structure of H+(S) is more complicated and still 
needs be explored. It is easy to show, however, that 
H+(S) contains Hi as a subset: Since Hi C H+(So) be- 
longs to the closure of D + (So) (and since So CS), it must 
also be included in the closure of D+(S). But since all 
points of Hi sits on CNGs, none of them belong to D+(S). 
Therefore all points of Hi must be located at the bound- 
ary of D_|_(S) (but not on S itself), namely on H+(S). 

The fact that Hi is contained in H + (S) guarantees the 
"causal protection" discussed in the previous sections. 
Consider a point P located in Hi but away from its in- 
tersection with H2 (i.e. at some < 9q — 9 c ). The bound- 
aries of D + (S) and D + (So) overlap in the neighborhood 
of N. Therefore, any sufficiently-small neighborhood of 
N in the closure of D + (S) is contained in the closure of 
D + (So), and is hence guaranteed to be regular. This en- 
sures that no singularity which evolves at the boundary 
of D + (S) may get close to P. Obviously this argument 
also applies to the entire CNGs located at 8 < 8q — 8 C . 



IV. THE INITIAL- VALUE SET-UP 

A. Basic strategy 

Our construction of the TM spacetime is formulated 
in terms of the corresponding initial data on the ini- 
tial hypersurface S. These initial data, which include the 
three-metric h a b and the extrinsic curvature K a b [25| , are 
constructed so as to satisfy the constraint equations (dis- 
cussed below). The evolution of geometry will in turn be 
determined by the evolution equations. The set of equa- 
tions relevant to our model is the Einstein-dust system, 
namely 

G af3 = 8nT at3 = 87reu Q ?/ , (8) 

where e is a scalar field and u a is a normalized vector 
field. These quantities correspond to the dust density 
and four-velocity, respectively. The Energy conditions 
are all satisfied if e > and u a is a timelike vector. The 
system (j8]), along with the initial data on S, uniquely de- 
termine the evolution of geometry (and matter) through- 
out D+(S) 1[H 

As was described in the previous section, the initial 
hypersurface S is composed of three parts: The inner 
core So, the envelope Si, and the external region S2. In 
So and S2 the initial data correspond to vacuum (i.e. 
e = 0). This [along with the mathematical properties of 
Eq. ©] guarantees that the evolving geometry will be 
vacuum throughout D + (So) and D + (S2). In the rest of 
D + (S), e will in general be positive [though it may also 
vanish in certain portions of J + (Si)]. 

The evolving vacuum metric is known explicitly 
throughout D+(So), it is given in Eq. ([2|) [or Eq. 
(HJ]. The vacuum solution in D+(S2), the external 
part of D+(S), is also known analytically: It is just 
the Schwarzschild solution. Our construction of S thus 
guarantees that the evolving spacetime in D+(S) will be 
asymptotically-flat and will admit future null infinity. In 
addition, the way we construct So (in particular the re- 
quirement 8q > 8 C ) guarantees that the conditions rel- 
evant to the central core region [e.g. features (5-7) in 
subsection II Al above] are satisfied. 

The "envelope" part of the evolving spacetime is the 
region between D+ (So) and D+(S2). It may be expressed 
as D + (S)nJ + (Si). The evolving metric in this part is 
not known to us. In particular, we do not know which 
kinds of singularities (if any) develop there, and where. 
This limits our present ability to analyze the full causal 
structure of our spacetime, e.g. whether an event horizon 
form, and where exactly is the CH (outside of D+(So)). It 
seems that a numerical solution of the evolution equation 
will be required in order to fill this gap. Nevertheless, 
the known analytic vacuum solutions throughout D + (So 
) and D + (S2), along with the properties of S and the 
initial data on it, guarantee that all the conditions (1- 
10) in subsection II Al are satisfied by our spacetime. 
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B. The constraint equations 

The initial data h a b and K a b, to be specified on S, 
are subject to four constraint equations, which corre- 
spond to four combinations of the Einstein tensor that 
are completely determined by h a b and K a b- Let N a be 
a normalized timelike vector (defined on S) orthogonal 
to S, and let x a be a set of three spacelike coordinates 
parametrizing S. Then the four constrained components 
of the Einstein tensor are 



G a = G aa N a — K a:l 



K, 



b : a 



and 



G = G a pN a N 13 = - \r<® + {K a a f - K ab K al 



where is the Ricci scalar associated with the 3-metric 
h a b, and a colon denotes covariant differentiation with 
respect to h a b- (Indices of K are rased and lowered with 
the three- metric h a b-) 

The Einstein equation G a p — &irT a p then imposes the 
three Momentum equations 



a:b 



K b . a = SwT aa N a = 8nT a 



and the Energy equation 

R {3) + {K*f - K ab K ah = \QnT af} N a N p = 16nf. 

In our dust model T a ^ = £u a u@ . To simplify the analysis 
we now choose the initial dust velocity u a to coincide 
with N a . Then f a vanishes (because N a = Q), and f = e. 
We obtain the momentum equation 



K b 



and the Energy equation 



K b 







K ab K ab 



(9) 



(10) 



In the internal part So the extrinsic curvature takes 
the simple form (|4l5p . We shall now adopt this same 
form of K a b for the envelope Si as well. Then in both So 
and Si the Momentum equation ^ is trivially satisfied, 
and the only nontrivial constraint equation is the Energy 
equation, which now reads 



(11) 



In the internal vacuum region So this equation reduces 
to 



2 



(12) 



and one can easily verify that the 3-metric Q indeed 
satisfies this relation. In the envelope region Si, we only 



require e to be non-negative, therefore the Energy equa- 
tion becomes an inequality 



(13) 



In the external vacuum region the extrinsic curvature 
will no longer take the form ^ . In this region the initial 
data will need to satisfy the Momentum equation ((9| as 
well as the vacuum Energy equation, namely 



i? (3) + (A^) 2 - K ab K ab = . 



(14) 



The initial data for the core So are given in Eqs. 
and ((4]). We shall now proceed to construct the initial 
data for Si and S2 as well. 



V. INITIAL DATA FOR THE ENVELOPE 

The envelope Si interpolates between the internal vac- 
uum core solution and the external Schwarzschild geom- 
etry. Throughout Si the extrinsic curvature is given by 
Eqs. (|4l5p . and the 3-metric only needs to satisfy the 
dust inequality (TIB"]) . 

This region will be further divided into three sub- 
regions: 

The inner part: A region which extends between the 
torus 9 = 9o (the inner boundary of Si) and a cer- 
tain larger torus 9 — 63 > 9 ; 

The outer part: An inhomogeneous spherically- 
symmetric dust solution near the outer boundary 
of the envelope. This sub-region extends between 
the 2-sphere R — R s bounding Si and a certain 
2-sphere R — R\ < R s (which contains the torus 
9 = 9 3 ); 

The intermediate part: A homogeneous dust solution 
which extends between the inner and outer parts 
(i.e. between the torus 9 = 63 and the 2-sphere 
R = Ri). 

Our construction of the initial data for Si will start from 
the inner part, proceed with the intermediate part, and 
conclude with the outer part. 



A. The inner part 

This region interpolates between the internal vacuum 
solution and a homogeneous dust solution. At its inner 
boundary the 3-metric is 



ds 2 = rl{d9 2 + sinh 2 9d<p 2 ) + L 2 a 



(15) 



and on approaching the outer boundary it will become 
the flat 3-metric (|2"4"|) below. A general form which covers 
both metrics (as well as the entire region in between) is 
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ds 2 = r 2 [d8 2 + F{8) 2 d^ 2 } + L(6, f) 2 d^ 2 



(16) 



In the various sub-layers of the envelope's inner part we 
shall make certain choices for the functions F(9, <p) and 
L(9, ip), so as to allow the transition from the metric (fl5|) 
to &ty without violating the positivity of the dust energy 
density. We require both F and L to be smooth (namely 
C' 00 - 1 ) and positive, which will ensure the regularity of 
the 3-metric (fT6|) . 

Our starting point is the 3-metric (| 1 5[) . namely 



L = const — Lq . F = sinh 8 , 



(17) 



at 8 < 8q. Then in the range 8q < 8 < 8\ (for some 
0i > 0'i nin > 0o, where 0™ in is given below), we still take 
L = Lq but allow a general function F(8) > 0, so the 
3-metric is 



ds 2 = rl\dd 2 + F{8) 2 dip 2 } + L 2 Q d(j) 2 . 
The Ricci scalar is then 

2 J> 
V F 



(18) 



= -2F.ee/irlF) 



Thus, apart from F(8) > 0, the function F(8) only needs 
to satisfy 

F,ee < F . 

(Note that this becomes an equality at the vacuum region 
< O , where F = sinh0.) 

We want the function F(8) to join smoothly on F = 
sinh0 at < O . We also want it, in view of Eq. (pM)) 
below, to join smoothly on F = 9 at 6 > Q\. It is not 
difficult to construct a function F{8) satisfying all these 
features — provided that 8\ is sufficiently large. In fact 
we take 



0i nin = 20 o + 1 . 

This construction of F(8) is described in Appendix A. 

Next, in the range 0i < < 02 (for any 02 > 0i), we 
take F(8) = 0, so the 3-metric is 



ds 2 = ii[dd 2 + 8 2 dip 2 ] + Lf,d(j) 2 



(19) 



This metric is flat, hence the dust inequality (| 1 31 ) is 
trivially satisfied. 

Next, in the range 02 < < 03 (for any 83 > 82} we 
take the 3-metric in the form 



ds 2 = r 2 [dd 2 + 8 2 d(p 2 } + L{8, <p) 2 dcj> 2 
One finds for this metric 

2 



(20) 



n -AL/L. 
V 1 



where A denotes the 2-dimensional flat-space Laplacian 
in polar coordinates, namely 

AL = L. ee + L,g/8 + L^/8 2 . 

The dust inequality (|13p then implies that L(8, tp) must 
satisfy 

AL < L . 

(along with L > 0). We now take L(8, if) in the form 
L(6,<p) =L + f(8)r 8 cos ip, 

where 

8 = (8-8 2 )/S8, 

and <50 = 03 — 02. The function /(0) may be any smooth 
function satisfying the following: 

(i) it joins smoothly on / = at < (corresponding 

to L(0,<p) = L Q at < 2 ), 

(ii) it joins smoothly on / = 1 at > 1 (corresponding 

to L(9, ip) = Lo + rod cos <p at > 03), 

(iii) it is monotonous in between. 
For this choice of L(0, <p) one finds 

V(0) , 0f"(8)\ 



AL = ro cos tp 



58 S8 2 



where a prime denotes here a differentiation with respect 
to 0. Using < 03 we obtain 



where and denote the maximal absolute values 
of df /dd and d 2 f/d8 2 , respectively, throughout the range 
< < 1. Noting that 



L(8, <p)>L - r 8 3 , 



(21) 



we conclude that the dust inequality AL < L is satisfied 
if Lq is greater than some minimal value 



- 7 '0#3 + [ 3 / r ' nax + 7^/, 



3 t n 

max 



(22) 



We shall thus take L to be > L™ m as required, hence 
the dust inequality is satisfied. This corresponds to 



I > ^min = 27TV / 3io 1 



(23) 



in terms of the parameters of the original core metric {T]) 
Note that Eq. (|2"Tj) and the above expression for L™ ln 
also guarantees that L(0, ip) is strictly positive through- 
out < 3 . 
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B. The intermediate homogeneous part 



Next, in the range 8 > 63 the 3-metric takes the form 

ds 2 = r 2 [<10 2 + 8 2 dtp 2 } + (L + r 8 cos tpfdcf) 2 . (24) 

(Recall that both ip and <j> admit a 2ir periodicity.) This is 
a flat metric in somewhat unusual coordinates. To bring 
it to a standard form we perform the following coordinate 
transformation: 



p = tqO cos (p + Lq , z = ro9 sin <p , 
and the metric becomes 



Note that 



ds 2 = dp 2 + p 2 d(f> 2 + dz 1 



P>L Q -r 6> L$ m - r Q t 



(25) 



hence p is positive on the torus 9 = 63 and in its neigh- 
borhood. 

In summary, outside the torus 9 = 63 the 3-metric is 
flat and the dust has a constant positive density, 



72irp 2 



The flat 3-metric and the uniform extrinsic curvature ((H) 
indicate that the initial data at 9 > 63 are just those 
of the (contracting) spatially-flat dust Robertson- Walker 
geometry. 



C. The outer part: inhomogeneous spherical dust 
geometry 

The outermost layer of the envelope Si will be con- 
structed to be a spherically-symmetric, inhomogeneous, 
dust region, which interpolates between the constant den- 
sity eo > at 9 > 9^ and the vanishing density at the 
Schwarzschild exterior. To this end we first transform 
the flat metric (|25p into standard spherical coordinates 
(R, 6, 4>) by 

z = R cos 8 , p — R sin <d . 
The 3-metric becomes 



ds 2 



dR 2 + R 2 dn 2 , 



(26) 



where dVi 2 is the unit 2-sphere, 

dVi 2 = dB 2 + sin 2 9# 2 . 

R may be expressed directly in terms of the original 
toroidal coordinates 9, tp: 



This implies an inequality 

L Q -r 9<R<L + r 9 . 



(27) 



We now truncate the flat metric (|26|) at a two-sphere 
R = Ri. We take 



i?i > L + 9 3 r 



(28) 



This ensures, by virtue of the second inequality in Eq. 
(f2"T)) . that the torus 9 = #3 is entirely contained at 
R < R\. The flat metric ([26]) thus holds throughout the 
region between the torus 9 — 6*3 and the sphere R = R\ 
surrounding it. 

The 3-metric at R > R± is assumed to be spherically 
symmetric, and we write it in the general form 



ds 2 = g RR (R)dR 2 + R 2 dfl 2 . 
It is convenient to substitute 

g RR = [1 - 2M(R)/R + R 2 /27p 2 Y 



The Ricci scalar is then 



2 



r co Z 4 dM 



4 dM 
R?~dR 



and Eq. (fTTjl yields 



1 dM 
4ttR 2 ~cIR 



(29) 



(30) 



(31) 



(32) 



Therefore, to ensure non-negative dust density, the func- 
tion M (i?) must be a monotonously-increasing (or at 
least non-decreasing) one. 

Consider next the boundary conditions on M(R). In 
the homogeneous dust region at R < R\ the 3-metric is 
flat, g RR — 1, which corresponds to 



M(R) = R 3 /54p 2 



(33) 



The external boundary of the spherically-symmetric dust 
region is at the two-sphere R = R2 for some R2 > Ri, 
where the dust solution is to be matched to a spherically- 
symmetric vacuum solution at R > R2 (see next section) . 
From Eq. (|32p this vacuum solution is characterized by 

M = const = m . 

Thus, M(R) is required to be a monotonously-increasing 
function which smoothly joins M = R 3 /54p, 2 at R < Ri 
and M — m at R > i?2. For later convenience we also 
demand 



M{R) < R 3 /Mp 2 



(34) 



It is straightforward to construct a function M(R) satis- 
fying all these requirements, for any given R2 > R\ and 
m > R\/5Ap 2 . We shall take 



R 2 



p 2 + z 2 = L 2 



(r #) 2 + 2L r 9 cos (p . 



R 2 >2m. 



(35) 



10 



From Eq. (M) it follows that 

1 - 2M(R)/R + R 2 /27fi 2 > 1 . 

This guarantees that the 3-metric (|29I30[) is regular 
throughout R 1 < R < R 2 . 

In the homogeneous dust region Eq. (I33|) implies 
2M/R = R 2 /27fj, 2 . Applying this to R = R 1} using 
Ri > 28 3 r [obtained from Eqs. (|22I28[) ] and 9 3 > 6 C > 1, 
one finds that 2M(fli) > R\. Namely, the dust solution 
includes spherical trapped surfaces at R = R\ and its 
neighborhood. 

Non-spherical modification 

The above construction of the inhomogeneous dust re- 
gion was spherically-symmetric. However it is easy to 
generalize it to obtain nonsphcrical configurations. This 
increases the space of solutions, and also allows for new 
kinds of causal structures. Since this modification goes 
beyond the main course of this paper, we shall only sketch 
it briefly here. 

In the first stage, one chooses the function M(R) such 
that at a certain region R a < R < Rb, for some i?i < 
R a < Rb < R 2 , it takes the form 

M{R) = aR 3 /5Afi 2 , (36) 

where < a < 1 is a fixed number. The 3-metric then 
becomes 

ds 2 = [l + ^^R 2 ]- 1 dR 2 +R 2 dn 2 . (37) 

This is a maximally-symmetric 3-metric of negative- 
curvature. Recalling the uniform extrinsic curvature ([4]), 
one realizes that the initial data at R a < R < Rb cor- 
respond to a contracting, unbounded (i.e. "k = — 1"), 
Robertson- Walker solution. Indeed from Eq. (f3"2"| the 
dust density is constant, e = aeo- One can easily ar- 
range the parameters a,R a ,Rb (and the mass function) 
such that the region R < R a is free of trapped surfaces 
(namely M(R) < R/2). 

Next, one picks a point P (corresponding to certain 
R,@,4>) on S somewhere at R a < R < Rb, and re- 
express the 3-metric in spherical coordinates centered 
on the point P (this is possible because the metric is 
maximally-symmetric). We denote these new spherical 
coordinates (R, 8, </>). The 3-metric in these new coordi- 
nates still takes the form (|3"7| , but with R,Q,tp replaced 
by R, 0, <f>, respectively. This 3-metric may be expressed 
by Eqs. (|29I30I36[) . with R and ft replaced by R and Cl, 
respectively. 

Finally one picks a two-sphere R = Rq around P which 
is entirely contained in R a < R < Rb. At R < Ro one 
modifies the mass function and picks a (monotonously 
increasing) smooth function M(R) at will. [Optionally 
one may also modify the extrinsic curvature, namely re- 
place Eq. Q by a more general spherically-symmetric 



form, in a manner described in the next section.] One 
then obtains an inhomogeneous dust solution at R < Rq. 

This modified solution at R < Ro is on itself spher- 
ically symmetric, but is not concentric with the spher- 
ical shells at e.g. R\ < R < R a . It therefore breaks 
the global spherical symmetry of the external parts of S, 
which would otherwise extend all the way from spacelike 
infinity to the trapped surfaces in the neighborhood of 
R = Ri . 

We mentioned above that this modification may allow 
for new types of causal structures. Here is one exam- 
ple: One can shape the inhomogeneous dust solution at 
R < Rq such that it will subsequently develop a naked 
shell focusing singularity (similar to those constructed in 
e.g. Ref. [l3} ). It is not difficult to arrange that this 
singularity will be globally naked. This spacetime will 
fail to be future asymptotically predictable 

VI. INITIAL DATA FOR THE EXTERNAL 
VACUUM REGION 

In the region R 2 < R < R3, for some R3 > R 2 , we set 
the 3-metric 

ds 2 = (1 - 2m/ R + R 2 / 27 (i 2 )- 1 dR 2 + R 2 dn 2 , (38) 

and the uniform extrinsic curvature ([4]). The Ricci scalar 
is R^ = — 2/9/i 2 and e vanishes. Note that the term 
1 — 2m/ R + R 2 /27n 2 was shown to be positive (if fact 
> 1) at R = i?2 and it is also monotonously increasing 
in R, hence it is positive throughout R 2 < R < -R3. 

Since the initial geometry in this range is both vacuum 
and spherically symmetric, it must correspond to that 
of a certain spherically-symmetric initial hypersurface in 
the Schwarzschild geometry. Although the Schwarzschild 
spacetime is asymptotically-flat, the initial 3-metric |38|) 
is obviously not. (In fact this three-metric is the same 
as that of a time-symmetric hypersurface in the de Sit- 
ter spacetime.) As it turns out, the initial data we have 
constructed in the range R 2 < R < R3 correspond to a 
" hyperbolic" (rather than time-symmetric) initial hyper- 
surface in the Schwarzschild geometry. Since we do want 
S to be asymptotically-flat, essentially what we need is 
to deform this initial hypersurface at R > R 3 (say), so 
that at large R it approaches a time-symmetric hyper- 
surface in Schwarzschild, namely one with vanishing ex- 
trinsic curvature and asymptotically-flat 3-metric. 

Here we describe the construction of S in terms of the 
initial data for h and K (rather than through its embed- 
ding in a given spacetime) . The discussion above makes it 
obvious, though, that in order to make S asymptotically 
flat we must relax the condition (01 on K a b at R > R3. 
(This would amount to " changing the embedding of S in 
spacetime" .) 

In the next two subsections we shall construct K a b and 
h a b, respectively, in the range R > R3. Our only pre- 
sumption is that both tensors are spherically-symmetric, 
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with h ab given by Eq. (|29p . The extrinsic curvature will 
be obtained from the momentum equation @, and the 
3- metric [namely the function g RR (R)] will in turn be 
derived from the vacuum Energy equation \l l 



A. The extrinsic curvature 

Being a spherically-symmetric tensor, we write K ab as 

K ab = K {R)h ab + 5K{R)n a n b , (39) 

where n a is the unit radial vector field. This expression 
must satisfy the momentum equation 



K" 



0. 



Owing to the linearity of this equation, we may consider 
the contribution of each term in Eq. (|39|) separately. 



Since K^h^ — 3Ko, the first term contributes —2Ko a . 
The contribution of the second term is 



[SK(R)n a n b ] , b - [SK(R)n b n b ] , a . 



(40) 



Since the gradient of SK is tangent to n, the contribution 
coming from the derivative of SK cancels out between the 
two terms in Eq. (|40"]) . Also (n b n b ) :a vanishes due to nor- 
malization. In addition the term n b n a . b vanishes by the 
geodesic equation, because n a is the tangent vector to a 
congruence of geodesies (the radial rays). The expression 
(|40|) therefore reduces to SKn a n b b , and the momentum 
equation becomes 



2K< 



0,a 



5Kn a n b b 



The angular components trivially satisfy this equation. 
In evaluating the radial component, a straightforward 
calculation yields 

n b b = (2/R)(g RR )~ 1 / 2 . 

Since n R = {gnu) 1 / 2 , the momentum equation reduces 
to the simple relation 



SK = R 



dK 
dR 



(41) 



Note that the choice ((4]) at R < R 3 corresponds to 
Kq(R) — const = k$. Then Kq(R) varies in the range 
R3 < R < R4, for some R 4 > R 3 . At R > R 4 , the 
outermost layer of S 2 , we choose K (R) = 0, which yields 
K a b = (this would correspond to a time-symmetric 
hypersurface in Schwarzschild) . 

Thus, in the transition region R 3 < R < R4 we take 
Kq(R) to be any smooth function which smoothly joins 
on K (R) = k at R < R 3 and on K (R) = at R > R 4 . 
The function SK(R) is then defined by Eq. (|4ip. hence 
the momentum equation is satisfied. 



B. The 3- metric 

The 3-metric at R > R3 will be determined from the 
vacuum Energy equation 

R^ + {K a a f-KtK b a = Q. 

With the substitutions (|39l4ip this equation becomes 



i? (3) = -2^0(3^0 + 2R 



dKo 
dR ' 

In the 3-metric (|29|) we now set 

g RR (R) = [1 - 2M(R)/R}~ 1 
The Ricci scalar is then found to be 

R 2 dR 



(42) 



(43) 



This, combined with Eq. (14"2")) . yields a closed expression 
for dM /dR in terms of Kq(R). After integration one 
obtains 

M(R) = m- R 3 (K ) 2 /2. 

The integration constant m is determined from the 
boundary condition at R = R 3 . Thus, in the range 
R3 < R < R4 (and, in fact, throughout R > R 2 ) the 
3-metric is (l29l) with 



g RR = [l-2m/R + R 2 K (Ry] 



21-1 



(44) 



Note that R 3 > R 2 > 2m, hence the 3-metric (|29|44l ) 
is regular throughout R > R 3 (its regularity at R < R 3 
was already established above). 

Finally, at R > R4 we have Kq = 0, hence the 3-metric 
and extrinsic curvature are 



and 



ds 2 = (1 - 2m,/Ry 1 dR 2 + R 2 dn 2 



K ab = 0. 



(45) 



(46) 



Obviously this corresponds to a time-symmetric initial 
hypersurface in a Schwarzschild spacetime with mass m. 



VII. DISCUSSION 

We have constructed an asymptotically flat spacetime 
which evolves from a regular partial Cauchy surface S and 
subsequently develops CCCs. The formation of CCCs 
takes place in a compact region of space; that is, S in- 
cludes a compact set So such that the CCCs form at the 
closure of D + (So). This central region is empty, and so 
is the external part of the spacetime. The intermedi- 
ate region (the "envelope") is made of dust with non- 
negative energy density. Although dust is not the most 
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realistic description of matter, it nevertheless provides a 
rather simple paradigm, which proved in the past to be 
useful in addressing various issues of principle in Gen- 
eral Relativity — e.g. gravitational collapse, formation of 
naked singularities, and cosmological models. Further- 
more, the system of dust+gravity is known to yield a 
well-posed initial- value problem [9(. 

The spacetime constructed in this way satisfies all the 
requirements (1-10) listed in subsection II At In particular 
it is smooth, asymptotically flat, and topologically trivial 
(to be precise, the initial hypersurface S is of topology 
R 3 ). It trivially satisfies the energy conditions (weak, 
strong, and dominant). 

The vacuum core metric was taken here to be the 
pseudo-Schwarzschild metric. We point out that we 
could also use the standard Schwarzschild metric (with 
the coordinate t identified on a circle), or even the 4- 
dimensional Misner space, for the core metric, and ob- 
tain a TM model with similar properties. However, it is 
only the pseudo-Schwarzschild metric which admits a ho- 
mogeneous initial hypersurface (namely r = 3(i) with a 
uniform extrinsic curvature ([4]). This simplifies the con- 
struction of the initial data for the envelope, because the 
momentum equation is automatically satisfied. With the 
alternative core metrics previously mentioned, the con- 
struction of initial data will be slightly more complicated. 

Several problems and important questions are still left 
open. Perhaps the most important one is the issue of 
stability. A stability analysis is beyond the scope of this 
paper, but there is comment worth noting. Although 
there are indications for classical and semiclassical insta- 
bilities in various TM models (see e.g. H^dHH), the 
robustness and effectiveness of these instabilities are still 
unclear [27] 29] [30] [32j] . Further research is required 
in order to assess the robustness and effectiveness of the 
various instability phenomena. The model constructed 
here may provide a more solid basis for a systematic sta- 
bility analysis. A few of its features that could be im- 
portant for a genuine stability analysis include: (i) hav- 
ing a regular initial hypersurface, (ii) asymptotic flatness, 
and (iii) admitting a well-posed system of evolution equa- 
tions. None of the previous TM models demonstrated all 
these properties. 

Two other important open questions should be men- 
tioned here: 

1. It may turn out that the evolving spacetime in- 
cludes a black hole, and all CCCs are imprisoned 
inside the event horizon. In such a case the for- 
mation of CCCs might still have crucial implica- 
tions to various aspects of the internal black-hole 
physics and geometry (e.g. singularity formation), 
but nevertheless the external universe will not be 
influenced. 

2. In our present construction, the initial data on S 
involve strong (though finite) gravitational fields. 
Is it possible to create a TM spacetime of this kind 
starting from weak-field initial data on some earlier 



initial hypersurface? Rephrasing this question: is it 
possible to create such a TM spacetime by sending 
weak gravitational waves (from past null infinity) 
and diluted dust shells (from past timelike infinity) 
in the inward direction? 

As a first step towards addressing these questions one 
may numerically evolve the initial data on S in both the 
future and past directions. Future time evolution will tell 
us whether a black hole forms, and if it does, whether 
it engulfs the CCCs. Past time evolution will proba- 
bly indicate one of the two possibilities: (a) the back- 
propagated fields disperses and weakens — the dust ex- 
pands and dilutes, and the gravitational field spreads to 
past null infinity as weak gravitational waves, or (b) the 
fields (back-) focus to form a white hole, a naked singu- 
larity, or pathologies of some other kind. Such numerical 
simulations will thus answer the questions (1,2) above at 
least with regards of the specific TM model constructed 
here. Note, however, that even if this specific model is 
indeed found to form a black hole in the future evolution, 
and/or a white hole in the past evolution, it leaves open 
the possibility that a modified TM model will be free of 
these undesired properties. I am not aware of any the- 
orem or argument which establishes a firm link between 
the formation of CCCs and the subsequent formation of 
a black hole, or the presence of a white-hole etc. to the 
past of S. 

One might hope to gain insight into these two questions 
by exploring the initial data on S for trapped and/or 
anti-trapped (i.e. "past-trapped") surfaces. Consider 
first the issue of trapped surfaces and their relation to 
black-hole formation in future time evolution. The ex- 
ternal spherically-symmetric vacuum region is free of 
trapped surfaces, because S2 is restricted to R > R2 and 
R2 > 2m. But as previously mentioned, the dust region 
in the neighborhood of R = Ri does include spherical 
trapped surfaces. However, the role of trapped surfaces 
as indicators for black hole formation is not so clear in 
our case. The theorems establishing the connection be- 
tween trapped surfaces and black hole formation assume 
either global hyperbolicity, or lack of CTCs, or asymp- 
totic predictability, or similar properties. Here none of 
these properties can be assumed apriori, especially be- 
cause the spacetime in consideration is guaranteed to 
develop CCCs and a CH. For example, on the basis of 
proposition 9.2.1 in Ref. Q, the occurrence of trapped 
surfaces on S basically tells us that one of the following 
two scenarios will be realized: (i) a black hole will form 
in D+(S) and engulf the CCCs, or (ii) the CH will extend 
to future null infinity, thus invalidating the condition of 
future asymptotic predictability. [331 ] 

The situation regarding anti-trapped surfaces seems to 
be different. If anti-trapped surfaces were found to be 
present on S, this would provide firm evidence that in 
past evolution, the field cannot just back-spread to infin- 
ity. It would indicate that D_(S) fails to be past asymp- 
totically simple, meaning that some pathology must have 
taken place in the past (prior to S): for example, a white 
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hole, or a naked singularity. Thus, the presence of anti- 
trapped surfaces on S would severely reduce the relevance 
of the spacetime in consideration as a physical model de- 
scribing the construction of a TM. (Presumably, the fu- 
ture "spacetime engineers" will not have white holes or 
naked singularities at their disposal.) 

Fortunately, in the specific model constructed here it 
appears that no anti-trapped surfaces exist on S. In re- 
gions So and Si, and also in the part R < R3 of S2, the 
simple form 0$ of the extrinsic curvature means that K a 
only has positive eigenvalues (triple fco > 0), which does 
not allow for anti-trapped surfaces. In the part R > R3 
of S2 the geometry is Schwarzschild with R > 2m, hence 
again there are no anti-trapped surfaces. 

Although the lack of anti-trapped surfaces is encour- 
aging, recall that it is a necessary but not a sufficient 
condition for the non-pathological asymptotic structure 
of D_(S). A numerical simulation of the initial data on S 
towards both the past and future directions could there- 
fore provide valuable insight into this issue of weak-field 
initial data prior to S, as well as into the problem of black 
hole formation in the future of S. 
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APPENDIX A 

We present here the construction of the function F(8) 
of Eq. (flgf in the range 6q < < 0i, where 6q and 0i 
are given numbers satisfying 80 > 8 C and 0j > 0™ 111 = 
28q + 1. Throughout this appendix a prime will denote 
differentiation with respect to 8. 

Recall the required properties of the function F(8): (i) 
It is strictly positive, (ii) F" < F, (iii) It is smooth in the 
range 0o < 9 < 0i, (iv) It joins smoothly on F(8) — sinh0 
at 9 < O , (v) It joins smoothly on F{8) = 6 at 6 > X . 

This construction is naturally divided into two stages: 
In the first one (stage A below) we construct a function 
F(9) which satisfies all the above requirements except 
that it is non-smooth (F' is discontinuous) at two points, 
denoted 80 and 8%, which satisfy 8q < 80 < 6\ < 6\. 
Nevertheless F(8) is continuous at 8 = 8q and 9 = 6\. 
Furthermore, we make sure that the "jump" in F'(8) is 
negative at both points #0 an d #1 (namely, the one-sided 
derivative in the direction 8 > #0,1 is smaller than the 
corresponding one in the direction 8 < #0,1 )■ Then in 
the next stage (referred to as stage B) the function F(8) 
is modified — smoothened — in narrow neighborhoods of 
both 8 = 8q and 9 = 6±, to obtain a smooth function 
F(8) in the entire range 8q < 8 < 8\ which satisfies all 
the above requirements (i-v). 



Here we shall describe in some detail the procedure 
comprising stage A. The main statement underlying 
stage B, namely, that a function with a "kink" (of the 
correct sign) can be smoothened without violating prop- 
erties (i,ii) above, is quite obvious, but the full presen- 
tation of the detailed smoothening procedure is lengthy. 
We shall therefore skip the detailed description of stage 
B 



Stage A: Constructing the rough function F(8) 



Let us define 



F(8) = a cosh(0-20 o ) 



(47) 



We take the function F(8) to be F(8) = F(8) at 8 < 
8 < 0i, along with F(8) = sinh0 at 8 < O and F(0) = 8 
at > 0i. The points 8q and 8\ will thus be the in- 
tersection points of F(8) with sinh 8 and 0, respectively. 
(The smooth matching at points 0o an d 0i is thus triv- 
ially satisfied, but as we mentioned above the challenge 
remains to arrange the smooth matching at 0o and 0i 
— the central task in stage B.) 
We take 

01 = 20o + 1 , 

hence 0o < 0i < 0i as desired. The parameter a will then 
be derived from the requirement of continuity at 8 = 8±, 
and subsequently continuity at 8 = 8q will determine the 
value of O . Note that Eq. g7]) satisfies F" = F (hence 
e vanishes — and condition (ii) above holds — throughout 
the range where F = F). 

Continuity at 8 — 81, namely F(8i) = 81, yields 



a = (20 o + l)/coshl . 



(48) 



Note that 



a > (20 c + l)/coshl ^ 2.35, 
hence in particular 

a > 1 . 



Also note that 



a < 20 o + 1 < e 2 



The parameter 0o is taken to be the point where F{8) 
intersects sinh 0, namely it satisfies 



a cosh(0o — 20o) = sinh 0o . 
This equation has a single real root: 



O = 0o + (1/2) In 



1 
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Note that 9q > 6q as desired. Later we shall also need 
the inequality 9 < 29$. To see this, one compares F(9) 
to sinh# at the two points 9 = 9q and 9 = 29q. In the 
former 

F(9 ) = a cosh 9 > cosh(9 > sinh0 o • 

On the other hand, at 9 — 29q one obtains 

F(26» ) = a = (26»o + 1 )/ cosh 1, 

which is to be compared to sinh(2#o)- One finds (e.g. 
numerically) that 

(26>o + l)/coshl < sinh(26> ) 

for any 9q > 9 C (the above inequality in fact holds for any 
9o greater than 0.57, whereas 9 C = 1.317). Since F(9) is 
smaller than sinh 9 at 9 = 29o but greater than sinh 9 at 
9 = 9a, the (single) intersection point 9q must be located 
in between, namely 

9 < 9 < 29q ■ 

Finally we compare the two one-sided values of F' at 
each of the matching points 9q and 6\. Starting at 9o, 
the directional derivative corresponding to 9 < 9q is 

F' = cosh <9 > , 

whereas the one corresponding to 9 > 9o is 

F' = F' — asinh((9 - 2 9 ) < , 



hence the jump in F' is negative. Consider next the two 
one-sided values of F' at 9\ . For the direction 9 > 9\ we 
have F' = 1, and for 9 < 9\ we have 

F' = F' = a sinh 1 . 

Substituting the value of a, Eq. (|4"5|) , we get (at 9 = 9\) 

F' = (29q + 1) tanh 1 > (29 c + 1) tanh 1 = 2.77 > 1 . 

We conclude that at both 9 = 9o and 9 = 9\ the jump 
in F' is negative (namely, the directional derivative cor- 
responding to 9 > #o,i is smaller than the corresponding 
one corresponding to 9 < #o,i), as desired. 

Stage B: Smoothening the rough function F(6) 

In the next stage one constructs the function F{9) in 
the range 9o < 9 < 9\ to be the same as F(9) except at 
two narrow ranges, one in the neighborhood of 9 = 9o and 
one in the neighborhood of 9 = 6\ , in which one replaces 
the non-smooth function F{9) by a smooth one. This 
must be done without violating the two equalities F > 
and F" < F. This procedure is quite straightforward 
though a bit tedious, and we shall skip the details here. 

It should be emphasized that this smoothening is only 
possible if the "jump" in F' is negative at both 9o and 
6\ (which was indeed shown above to be the case). For, 
only in this case F{9) satisfies the condition (ii) above in 
a distributional sense. 
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